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We compute the DC conductance for a homogeneous sine- 
Gordon model and an impurity system, which in the con- 
formal limit can be reduced to a Luttinger liquid, by means 
of the thermodynamic Bethe ansatz and standard potential 
scattering theory. We demonstrate that unstable particles 
and resonances in impurity systems lead to a sharp increase 
of the conductance as a function of the temperature, which is 
characterized by the Breit-Wigner formula. 
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I. INTRODUCTION 

In the context of integrable quantum field theories in 
I+I space-time dimensions a large arsenal of extremely 
powerful non-perturbative techniques has been developed 
over the last two and a half decades. The original moti- 
vation to treat these theories as a testing ground for re- 
alistic theories in higher dimensions is nowadays supple- 
mented by the possibility of direct applications, since the 
nanotechnology has advanced to such a degree, that one 
dimensional materials, i.e. quantum wires, may be real- 
ized experimentally. A quantity which can be measured 
directly [1] is the conductance through the quantum wire. 
There exist also already various proposals [2,3] of how 
to obtain this quantity from general non-perturbative 
techniques, such as the thermodynamic Bethe ansatz 
(TBA) [4,5] and the form factor approach [6] to com- 
pute the current-current two-point correlation functions 
in the Kubo formula [7] . Here we want to concentrate on 
the former approach. Whereas in [2] the emphasis was 
put on reproducing features of quantum Hall systems and 
the authors appealed extensively to massless models, we 
want to treat here in contrast systems which are purely 
massive. In particular we want to investigate how the 
properties of unstable particles and impurity resonances 
are reflected in a possible conductance measurement. 

II. FROM CONDUCTANCE TO MASSES OF 
UNSTABLE PARTICLES 

The direct current / through a quantum wire can be 
computed simply by determining the difference of the 
static charge distributions at the right and left constric- 
tion of the wire, i.e. / = Qr — Ql- This is based on 
the Landauer transport theory, i.e., on the assumption 



[2,3], that Q(t) ~ {Q R - Q L )t — {p R - p L )t, where the 
ps are the corresponding density distribution functions. 
For more details and a comparison with the Kubo for- 
mula see [8]. Placing an impurity in the middle of the 
wire, we have to quantify the overall balance of particles 
of type i and anti-particles i carrying opposite charges 
qi = —qi at the end of the wire at different potentials. 
This is achieved once we know the density distribution 
pi (9,r,[ii) as a function of the rapidity 9, the inverse 
temperature r and the chemical potential In the de- 
scribed set up half of the particles of one type are already 
at the same potential at one of the ends of the wire and 
the probability for them to reach the other is determined 
by the transmission amplitude |Tj (9)\ through the impu- 
rity. Therefore 

I = J2li(r^i) (1) 

i 

= E /^f U ( 6 > r >rf) - Pi vW)) \T? [0)\] ■ 

i 

By definition the DC conductance results as 

G M = E, = E, ^ w) M (2) 

and is of course a property of the material itself and a 
function of the temperature. In general the expressions 
in (1) tend to zero for vanishing chemical potential such 
that the limit in (2) is non-trivial. 

Let us now compute the density distribution by means 
of the thermodynamic Bethe ansatz. As was pointed out 
in [9], the TBA-equations for a bulk system and a sys- 
tem with a purely transmitting defect are identical. This 
is due to the fact that in the thermodynamic limit the 
number of defects is kept fixed and is therefore insignif- 
icant in thermodynamic considerations. Therefore the 
same equations also hold when we allow the impurity to 
be such that transmission and reflection are simultane- 
ously possible. We recall the main equations of the TBA 
analysis which are directly relevant in this context, see [4] 
for more details and in particular for the introduction of 
the chemical potential see [5]. For a detailed derivation 
of the TBA equations in this context see [8] . The main 
input into the entire analysis is the dynamical interaction 
encoded into the scattering matrix Sij (9) of two particles 
of masses mi and rrij and the assumption on the statis- 
tical interaction which we take to be fcrmionic. As usual 
[4,5], by taking the logarithmic derivative of the Bethe 
ansatz equation and relating the density of states Pi{9) 



1 



for particles of type i to the density of occupied states 
Pi(0) one obtains 



PiiB^in) = ^ coahe + Hfrv * ■ ^ 

By (f*g)(9):= I / (2ir) f d6' f (9 - 9')g(9') we de- 
note the convolution of two functions and ifij (9) = 
—id In Sij (9) jd9. The mutual ratio of the densities serves 
as the definition of the so-called pseudo-energies e$ (9) 



-£i(6,r,Hi) 



P T i{Q,r, fij) = 
Pi(6,r,Hi) l + e - e -( .^)' 



(4) 



which have to be positive and real. At thermodynamic 
equilibrium one obtains then the TBA-equations, which 
read in these variables 

rrrii cosh 9 = £j(0,r,/Xj) + rpi + ^2 [<Pij * Lj](6) , (5) 

where r = m/T, mi — » mi/m, pi — > Pi/m, Li(6,r,pt) — 
ln(l + e _£i ( e,r ^ i )), with m being the mass of the lightest 
particle in the model and T the temperature. It is im- 
portant to note that pi is restricted to be smaller than 
1. This follows immediately from (5) by recalling that 
£i > and that for r large Si(9,r,pi) tends to infinity. 
As pointed out already in [4], here just with the small 
modification of a chemical potential, the comparison be- 
tween (3) and (5) leads to the useful relation 



Pi(0,r,fj,i) = — 



1 / dsi(9,r,fii) 



dr 



+ P, 



(6) 



The main task is therefore to solve (5) for the pseudo- 
energies from which then all densities can be recon- 
structed. In general, due to the non-linear nature of the 
TBA-equation, this is done numerically. However, in the 
large temperature regime one may carry out various an- 
alytical approximations. For large rapidities and small r, 
one [4] can approximate the density of states by 

ft (g,r, W )~-ell~— e(6) - , (7) 

where e(0) is the step function. To obtain this we as- 
sume in (4) that in the large rapidity regime p r i {9 1 r 1 pi) 
is dominated by (7) and in the small rapidity regime by 
the Fermi distribution function, therefore 

pK9,r,p l )^^-e(9)-^\n[l + cxp(-e t (9,r,p l ))} . (8) 

Using this expression, we approximate the current in (1) 
and for pf- = —pf = V/2 the conductance results to 



lim G t (r ) ~ 
r^O w 2nr 



d9 



l + exp[ £i (0,r,O)] 



de t (9,r,V/2) 
dV 



d[<9)\Tm? 



(9) 



In order to evaluate (1) and (9) it remains to specify 
how to compute the transmission amplitude. In principle 
this can be done by exploiting the factorization equations 
which result as a consequence of integrability. However, 
for systems with a diagonal bulk S-matrix these equations 
are not restrictive enough and we will below simply use 
a free field expansion and proceed in analogy to stan- 
dard quantum mechanical potential scattering. Having 
obtained Tj(9) and Rj{9) one can construct the equiva- 
lent quantities for multiple defects from these functions. 
Here we are particularly interested in a double defect. 
Placing the two defects of the same type at x\ = 0, X2 = y 
the total transmission amplitude Tj can be build up from 
the ones of a single defect as [10] 



Tj(9) 



Tf(9) 



1 - R 2 A9) cxp(i2ysinht9) 



(10) 



Having assembled all the ingredients for the computa- 
tion of G we turn to the question of how the properties of 
unstable particles are reflected in this quantity? Assum- 
ing that (9) possesses a resonance pole at 9r = cr — ia, 
the Breit-Wigner formula [11] allows to determine the 
mass Mg and the decay width Tg of an unstable particle 
of type c 



2Mf = vV + 7 2 + 7 > 2(m 4 + m 3 ) 5 



(11) 



r~/2 = \/-f 2 + j 2 - 7 > 4miTOj(l - cosher cosct), (12) 

where 7 = mf + m 2 + 2mimj cosh a cos a and 7 = 
2rriimj sinh \a\ sina . The thresholds in (11) and (12) re- 
sult from energetic reasons [12]. We will now demon- 
strate that besides unstable particles also resonances 
in impurity systems can be described by means of the 
Breit-Wigner formula. An important consequence of 
(11) is that we can approximate the mass in there by 
M| w l/2mjrnj(l+cos<7) exp \a\ for large a. Then under 
a renormalization group flow Mg — ► rc-Mg the quantity 



Mg 



r^e 172 ^ 2 remains invariant. Once the 



unstable particle can be created, it can participate in the 
overall conductance and one should observe an increase 
at T c in G related to this process. For this interpre- 
tation to hold, we should observe the following scaling 
behaviour of the conductance 



G(r 1 c ,a 1 )=G(r 2 c ,a 2 ) for r*W 2 - r 2 c e 



01 li 



(13) 



Surely rc is not sharply defined, but taking for instance 
the middle between the beginning and the end of the 
onset seems reasonably well identifiable. Here rc is the 
inverse of the critical temperature rc = rn/Tc at which 
the unstable particle for fixed a is formed. This means 
the identification of the onset in a conductance measure- 
ment will provide rc, such that for given a the mass of 
the unstable particle can be deduced. 
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III. HOMOGENEOUS SINE-GORDON MODEL 



IV. FREE FERMION WITH IMPURITIES 



The £[/(3)2-homogeneous sine-Gordon (HSG) model 
is the simplest of its kind and contains only two self- 
conjugate solitons, which we denote by "+", "— ", and 
one unstable particle, which we call c. The correspond- 
ing scattering matrix was found [13] to be S±± = — 1, 
S±^(9) = ^ tanh (9 ± er — in/2) /2, which means the res- 
onance pole is situated at 9r = =f<t — i7r/2. Stable bound 
states may not be formed. It is known [10], that inte- 
grable parity invariant impurity systems with a diagonal 
bulk S-matrix, apart from S = ±1, do not allow simul- 
taneously non-trivial reflection and transmission ampli- 
tudes. This statement can be extended to the parity 
violating case [15]. We treat therefore (1) for a trans- 
parent defect, i.e. |T| = 1. The results for the conduc- 
tance after solving numerically the TBA equations for 
Hn , = —fiL = 0.25 are depicted in figure 1. 




log(r/2) 



FIG 1: Conductance G for the S'(7(3)2-HSG-model as a func- 
tion of log(r/2) = log(™/2T) for various values of the reso- 
nance parameter a. 

To carry out the limit /i — > is rather complicated 
when one does not have an explicit analytic expression 
at hand as in our case. However, we can take the result 
for finite a as a very good approximation, since we ob- 
serve that G(r)//j, ~ const for small r. We observe the 
onset of the unstable particle in form of a relatively sharp 
increase in G and in particular the validity of (13). The 
interpretation is clear: Only when we reach an energy 
scale at which the unstable particle can be formed it can 
participate in the conducting process. All this informa- 
tion is encoded in the density (9, r,/ij). Also the bound 
in (11) is respected. Computing now Ei(9, 0, 0) in a stan- 
dard TBA fashion, e.g., [16], we predict analytically the 
plateaux from (9) at 2(1 + v / 5)/(5 + y/E)n and at 1/27T. 
The latter value is obtained from the fact that in the re- 
gion in which a >• — 2 log(r/2), the system can be viewed 
as consisting out of two free Fermions such that (9) gives 
the quoted value. 



The continuous version of the 1+1 dimensional Ising 
model with a line of defect was first treated in [17]. 
Thereafter it has also been considered in [18,10] and [19] 
from a different point of view. In [17,18,10] the impurity 
was taken to be of the form of the energy operator and in 
[19] also a perturbation in form of a single Fermion has 
been considered. Here we also include a further type of 
defect. 

Let us consider the Lagrangian density for a complex 
free Fermion tp with £ defects 

C = <Ki7% - m)4, + J"" S(x - x n )V n $, i/)) , (14) 

* — ^n— 1 

where we describe the defect by the functions D„(^>, tp), 
which we assume to be linear in the Fermi fields. In 
the following we will restrict ourselves to the case 1 = 2 
with x n — ny and T> n (ip,ip) = T>(^),ip). We compute 
the transmission amplitude as indicated in [18,10,19], 
namely by decomposing the solution to these equations 
as ip(x) — 8(x) ip+(x) + Q(—x) ip-(x) and substituting 
them into the equations of motion. This way we obtain 
the constraints 



(15) 



Using now the standard Fourier expansion for a com- 
plex free Fermi field, the transmission and reflection am- 
plitudes can be read off componentwise from (15) as 
the coefficients of a]_{9) = R 3 {9)a] _{-&), a]_{9) = 

Tj{6)a] A6), etc. 

Recalling now that for the free Fermion the TBA- 
equations are simply solved by Ei(9, r, /ij) = rrrii cosh# — 
r/ij , we compute 



(16 



cosh 9 \T t {8)\ 2 
1 + cosh(r cosh t 



(16) 



To proceed further we have to specify the impurity. 



A. The energy operator defect, T>(ip, if)) = gipip 



From (15) we compute 
R j (6,B) = R J (9,B) 
T j (9,B)=T J (9,B): 



i sin B cosh 9 
sinh 9 + i sin B ' 
cos B sinh 9 
sinh 9 + i sin B ' 



(17) 
(18) 



where we used a common parameterization in this con- 
text sini? = — Ag/(A + g 2 ). The expressions Rj(9, B) and 
Tj(9, B) coincide with the solutions found in [10], which, 
however, in general does not correspond to taking our 
particles simply to be self-conjugate, since we use Dirac 
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Fermions. Using (17) and (18) we compute with (10) the 
conductance for a double defect with varying distance y. 
The results of our numerical computations for the con- 
ductance are depicted in figure 2. 

In the high temperature regime we can confirm once 
more these data by some analytical computations. From 
(9) we obtain 



1 / cos 2 B 
2vr 1 1 + sin 2 B 



<G(r,y) < 



1 

2^ 



(19) 



The lower bound becomes exact for y/r < 1. For B 
0.51 the values 0.0602 are well reproduced in figure 2. 
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FIG 2: Conductance G for the free Fermion with double en- 
ergy defect at distances y as a function of the inverse temper- 
ature r = m/T. 

We observe a similar type of behaviour as in the pre- 
ceding section and denote again the point of onset in the 
conductance by tq- Then, we deduce from our data that 
the following scaling relations 



G(r 1 c ,yi) = G(r 2 c ,y 2 ) for 



rcVi 



r cV2 



(20) 



Comparison with (13) suggests that we can relate the 
distance between the two defects to the resonance pa- 
rameter as a = 2 ln(const/y). However, despite the fact 
that the net result with regard to the conductance is the 
same, the origin of the onset is different. Whereas in the 
previous section it resulted from a change in the density 
distribution function it is now triggered by the structure 
of |T(0)|. Since p r keeps its overall shape and just moves 
its peak with varying temperature, the onset has to oc- 
cur when |T(0)| reaches its maxima. Using (17), (18) and 
(10), it is easy to verify that \f{6 = \n[{2n+l)n / y})\ « 1. 
Drawing now an analogy to the scattering matrix, this 
value plays the same role as Or and we therefore identify 



ln[(2n + l)%/y] ■ 



(21) 



Having fixed the resonance parameter a we may, in view 
of (20), relate the temperature to the mass scale of the 



unstable particle, associated now to the resonance, analo- 
gously as in the discussion after (13). However, there are 
some differences. Whereas in the HSG-model the onset 
is attributed to a single particle, the effect for the dou- 
ble defect system is attributed to several resonances. We 
identify a rs <tq +<ji • The other difference is that y is now 
a measurable quantity, such that a in (21) can be exper- 
imentally determined. On the other hand the sigma in 
(13) is usually a free parameter in the HSG-type models. 
Let us now verify our observations for a different type of 
defect. 



B. Luttinger type liquid, T>(tp,ip) = tp(gi + g2 r y°)tp 

There exist various ways to realize Luttinger type 
liquids [20]. Taking the conformal limit of the defect 
T>("tp,tjj) = ijj(gi + (727°)^: we obtain an impurity which 
played a role in this context [21] when setting the bosonic 
number counting in there to be one. Analogously to the 
previous sections we compute the related transmission 
and reflection amplitudes 



Tj(0,gi,g 2 ) 



4? (52 - g\ cosh( 



4«(.9i - 32 cosh6>) + (4 + g\ - g\)mti\d 

(4 + ffi-ff 2 )sinhfl 

4i(<jri - g 2 cosh6>) + (4 + g\ - g%) sinhfl 



The expressions for the particle j are obtained by replac- 
ing g\ — > —g\. The results of our numerical computation 
for gi — 0.7 and 172 = 0.2 depicted in figure 3 confirm the 
same physical picture as outlined in the previous sub- 
section. Our analytical prediction for the lowest plateau 
from (9) is 0.0324. 




FIG 3: Conductance G for the free Fermion with two defects 
D„(i/i,i/)) = jp(gi + 327°)^ at distance y as a function of the 
inverse temperature m/T. 
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V. CONCLUSIONS 



By using the TBA to compute the density distribution 
function and relativistic potential scattering theory to 
determine the transmission amplitude, we evaluated the 
DC conductance by means of equation (1). We demon- 
strated that the sharp increase of the conductance as a 
function of the temperature can be attributed to the pres- 
ence of unstable particles in the HSG models or likewise 
to a resonance of a double defect system. 
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